Halos and voids in a multifractal model of cosmic structure 
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ABSTRACT 

On the one hand, the large scale structure of matter is arguably scale invariant, and, on the 
other hand, halos and voids are recognized as prominent features of that structure. To unify both 
approaches, we propose to model the dark matter distribution as a set of fractal distributions of 
halos of different kinds. This model relies on the concept of multifractal as the most general scaling 
distribution and on a plausible notion of halo as a singular mass concentration in a multifractal. 
Voids arise as complementary to halos, namely, as formed by regular mass depletions. To provide 
halos with definite size and masses, we coarse-grain the dark matter distribution using a natural 
length derived from the lower scaling limit. This allows us to relate the halo mass function to the 
multifractal spectrum. Hence, we find that a log-normal model of the mass distribution nicely fits 
in this picture and, moreover, the Press-Schechter mass function can be recovered as a bifractal 
limit. 

To support our model of fractal distributions of halos, we perform a numerical study of the 
distribution produced in cosmological iV-body simulations. In the Virgo A-CDM GIF2 simulation, 
we indeed find fractal distributions of halos with various dimensions and a halo mass function of 
bifractal type. However, this mass function is just beyond the Press-Schechter's range, and we 
interpret it instead as caused by the undersampling of the distribution at the scale of halos, due 
to discretization. 



Subject headings: 
ods: statistical 

Introduction 



cosmology: large-scale structure of the universe - galaxies: clusters: general - meth- 



Current research on the large scale structure 
of matter is concerned with two fundamental con- 
cepts: halo and void. The concept of halo applies 
to the dark matter distribution and is meant to 
define the basic unit of structure. Therefore, dark 
matter halos have theoretical rather than observa- 
tional basis. In particular, the concept of halo is 
useful in the analysis of TV-body cosmological sim- 
ulations. In contrast, the concept of void arose in 
observations of the galaxy distribution, although, 
recently, it has also been applied to the dark mat- 
ter distribution and gained some theoretical basis. 

On a more fundamental level, the principle of 
scale invariance has been often applied to the 
study of the large scale structure of matter. Scale 
invariance is a general symmetry, with applica- 



tions in various branches of physics and other sci- 
ences. Its application in cosmology is old, dating 
back to hierarchical models of the Universe de- 
vised before the discovery of General Relativity 
(Mandelbrot, 1977). It became part of standard 
cosmology when the analysis of galaxy catalogues 
proved that the two-point correlation function is 
a power law on scales of several Mpc. A recent 
review of various ideas and models in cosmology 
based on scaling laws is given by Jones et al (2004) . 
Of particular importance are fractal models (Sy- 
los Labini, Montuori & Pietronero, 1998). Mul- 
tifractal models (Parisi & Frisch, 1985; Halsey et 
al, 1986) are an improvement of fractal models 
and were introduced in cosmology by Pietronero 
(1987), Jones et al (1988) and BaUan & Scha- 
effer (1988). Multifractal analyses of large-scale 
structure were initially applied to the distribu- 



tion of galaxies (Dominguez-Tenreiro & Martinez, 
1989; Martinez, Jones, Dominguez-Tenreiro & van 
de Weygaert, 1990) and later applied to TV-body 
simulations (Colombi, Bouchet & Schaeffer, 1992; 
Valdarnini, Borgani & Provenzale, 1992; Yepes, 
Dominguez-Tenreiro & Couchman, 1992). 

We shall take scale invariance of the dark mat- 
ter distribution in the nonlinear regime as our 
working hypothesis. This hypothesis leads es- 
sentially to a multifractal model (with transition 
to homogeneity). Since the linear regime de- 
parts from the initial scale invariance (Harrison- 
Zeldovich spectrum), the initial conditions for 
nonlinear evolution arc not scale invariant. How- 
ever, the dynamical equations for cold dark mat- 
ter are scale invariant. Therefore, it is reasonable 
to assume that the dynamics is driven to a scale- 
invariant attractor that is independent of the ini- 
tial conditions, as is normal in nonlinear systems. 

The power-law form of the galaxy-galaxy cor- 
relations can be taken as evidence of scale invari- 
ance, although there may be deviations from pure 
power laws. In particular, Zehavi et al (2004) 
find small but systematic deviations in the pro- 
jected correlation function of a sample from the 
Sloan Digital Sky Survey. They relate those devi- 
ations to the deviations from a power-law correla- 
tion function found by Jenkins et al (1998) in A- 
CDM simulations. Hence, they interpret them as 
caused by the transition from a two-halo regime on 
large scales to a one-halo regime on small scales. 
We shall discuss an alternative interpretation at 
the end of this paper, after we have introduced 
and explained the multifractal model of halos. Our 
interpretation is fully consistent with scale invari- 
ance. 

Halo models of large-scale structure are now 
wen developed (Cooray & Sheth, 2002). They 
assume that the full dark matter distribution is 
composed of two parts: the dark matter distri- 
bution within halos and the distribution of halos 
themselves. The matter distribution within halos 
is considered to be in virial equilibrium. This inner 
distribution has received much attention. Usually, 
it is assumed to be relatively smooth, and it is 
useful to define the average "profile" of a single 
halo. The first approximation is a power-law pro- 
file (Peebles, 1974; Sheth & Jain, 1997; Murante et 
al, 1997). However, more refined profile models in- 
volve a crossover between two power laws (Cooray 



& Sheth, 2002). 

Regarding the distribution of halos, halo cen- 
ters have been distributed randomly in some mod- 
els (Peebles, 1974; McClelland & Silk, 1977; Sheth 
& Jain, 1997; Murante et al, 1997) but clustering 
in the distribution of halos has also been consid- 
ered (Mo & White, 1996; Sheth & Tormen, 1999). 
Clustering of halos is to be expected and, in par- 
ticular, the clustering of halos can be of fractal 
type. In fact, a power-law halo profile and a frac- 
tal distribution of halos arc just two aspects of a 
new formulation of the multifractal model (Gaite, 
2005- A). In this formulation, a general multifrac- 
tal is considered as a set of fractal distributions of 
halos, such that similar halos have a fractal distri- 
bution with a given dimension. 

However, the separation of the full matter dis- 
tribution between a distribution within halos and 
a distribution of halos is not meaningful if the full 
distribution is scale invariant. A power-law profile 
actually prevents one from assigning definite sizes 
or masses to halos, which must be associated, in 
principle, with point-like mass singularities (Gaite, 
2005- A). But the dark matter distribution is not 
scale invariant on ever decreasing scales. Thus, 
the breaking of scale invariance induces a change 
in the nature of the matter distribution that allows 
us to assign halos a definite size. 

More in general, we will see that the fine struc- 
ture intrinsic to mathematical multifractals on 
ever decreasing scales gives rise to some counter- 
intuitive and even paradoxical geometrical no- 
tions. Some of these notions are relevant to de- 
scribe the dark matter distribution, in spite of its 
not being scale invariant on ever decreasing scales. 

Thus, our definition of halo differs from the 
usual definition: we understand that a natural def- 
inition of halo in a multifractal model with a nat- 
ural coarse-graining scale assigns halos equal size, 
namely, the size of that scale. A similar definition, 
namely, based on a uniform size, was proposed by 
Vergassola et al (1994) to define a mass function 
for the structure produced in the adhesion model. 
In contrast, the most popular definition of halo 
is based on the spherical collapse formalism (Mo 
& White, 1996; Sheth & Jain, 1997, Cooray & 
Sheth, 2002), and it assumes that halos are given 
by a density-contrast threshold {Sp/p — 178). 
This definition of halo constitutes the basis for 
the halo finders that use the "friends-of-friends" 



(FoF) technique (Davis et al, 1985). Nevertheless, 
Valageas, Lacey & SchaefFer (2000) have shown 
that both choices, namely, halos given by a den- 
sity condition or equal-size halos, can be consid- 
ered adequate from a theoretical point of view. 
Our multifractal model of halos agrees with their 
conclusion, for its crucial feature is the pattern in 
which the mass concentrates rather than the divi- 
sion of mass among individual halos, which can be 
made in various manners. Of course, the simplest 
manner is to take equal-size halos. 

The notion of void as basic element of the large 
scale structure is now well established, but its defi- 
nition is more imprecise than the definition of halo. 
The original definition of void as a totally empty 
region, suitable for the galaxy distribution, is now 
being replaced with definitions that allow for par- 
tial emptiness, after having been admitted that 
dark matter and even low luminosity galaxies are 
present in the initial voids. In fact, recently, Shan- 
darin, Shcth & Sahni (2004) have defined void as 
an underdense connected region, such that voids 
are complementary to superclusters. We propose 
here a definition that is related with the definition 
of Shandarin, Shcth & Sahni (2004), but which 
hinges on the nonlinear multifractal nature of the 
distribution, unlike theirs. 

We begin with a brief review of of multifractals, 
focusing on the concept of mass concentrations 
(singularities) and the realization of scale invari- 
ance. We introduce the key concept of multifrac- 
tal spectrum and, in particular, the bifractal spec- 
trum, as well as the method of calculation of the 
spectrum from correlation moments. Halos can be 
identified with singular mass concentrations, with 
power-law profile, whereas voids can be identified 
with regular mass depletions. Useful computa- 
tional notions of halo and void result when we 
consider a (natural) coarse-graining scale. This 
leads us to a general study of the coarse-grained 
probability distribution function and its connec- 
tion with the lognormal model. Further connec- 
tion with the Press-Schcchter mass function is also 
possible. Then we study what implications a mul- 
tifractal model has for the distribution of halos (as 
we have defined them). Finally, we test all these 
notions against cosmological iV-body simulations. 



2. Multifractals and mass concentrations 

Multifractal measures appear frequently in 
physics. In particular, they appear frequently 
as attractors of dynamical systems (Halsey et al, 
1986; Falconer, 2003). In dynamical systems, the 
natural measure of some region represents the por- 
tion of time that the system spends in it. For many 
systems, one indeed finds that, at long times, the 
trajectory is confined to a small region of the 
available space (an attractor), but the time that 
it spends in any subregion of it is very variable. 
In other words, the measure is very irregular and 
this non-uniformity encodes further structure that 
defines the attractor. 

The formation of large scale structure is given 
by the cosmological Newton-Poisson equations 
(Peebles, 1980), which constitute a dynamical sys- 
tem. Therefore, to justify the relevance of mul- 
tifractal analysis in the description of large scale 
structure, we may adopt the natural view that 
this structure is just an attractor of that many- 
dimensional dynamical system. Or we may argue 
that we expect a scale-invariant structure, but in 
which the realization of scale invariance is non- 
trivial, as in general multifractals. At any rate, 
the final justification must rest on a successful 
comparison with observational data or results of 
cosmological simulations. 

As said above, multifractal measures represent 
mass distributions spread according to highly ir- 
regular patterns, that is, with mass concentrations 
of very different magnitude (Falconer, 2003). This 
magnitude is referred to as "strength" and defined 
by the local dimension a (also called Lipschitz- 
Holder exponent in the mathematical literature); 
namely, 



a{x) = lim 



log?7i[i3(a;,r)] 
logr 



(1) 



where m[B{x, r)] is the mass in the ball of radius r 
centered on the point x. We also write, informally, 
m[B{x,r)] ^ 7."(^). \Ye only consider points x 
such that 'm[B{x, r)] does not vanish for any r > 0, 
to prevent the divergence of a{x). These points 
form the support of the measure. 

In a regular mass distribution, a ~ 3 (con- 
stant). So a multifractal is a singular mass distri- 
bution and its mass concentrations are also called 
singularities (of various strengths) . On the other 



hand, an ordinary self-similar fractal can be con- 
sidered endowed with a uniform mass distribution 
over it, such that a < 3 is the constant fractal 
(Hausdorff-Besicovitch) dimension. Thus, in the 
context of multifractals, ordinary self-similar frac- 
tals are called monofractals (or unifractals) . Full- 
fledged multifractals possess a range of exponents 
a (or strenghts), namely, < ttmin < o; < amax- 
Every set of points in which a takes a definite 
value is a fractal set (in the sense of having non- 
trivial fractal dimension). Therefore, a multifrac- 
tal can be considered as a set of fractals containing 
the various mass concentrations. The multifractal 
spectrum f{a) is the function that gives the fractal 
dimension of the set of points with exponent a. 

Regarding multifractal spectra, we should men- 
tion that there are cases in which a is not strictly 
constant but which cannot be considered full- 
fledged multifractals. An example is a finite 
number of (isolated) concentrations of different 
strength. This is a regular distribution except for 
the finite number of singularities. It has interest in 
cosmology, as the basis of the model of randomly- 
placed power-law halos (Gaite, 2005- A). The lo- 
cations of these halos may also be correlated, so 
that the model actually becomes related to the 
halo model proposed here, as we shall see. 

2.1. Correlation moments and multifractal 
spectrum 

To analyze the structure of multifractals, it is 
useful to introduce the scale-dependent correlation 
moments 



M,{r) 



jdH^)m[BMr\ 



(2) 



Ml is the total mass, which we normalize to 
one. With this normalization, a multifractal mea- 
sure can be interpreted as a probability measure. 
M2{r) is the usual mass-radius relation or two- 
point correlation integral (the integral of the two- 
point correlation function). Larger integer val- 
ues q — 3, 4, . . . are also employed in the analy- 
sis of the large scale matter distribution, but they 
are increasingly difficult to estimate. In addition, 
since multifractals are singular non-uniform distri- 
butions, integer values of q are not sufficient and 
we have to consider the full set of moments Mq (r) 
for — oo < g < oo. We can then define the function 
that gives the scaling behaviour of this full set of 



moments, namely. 



lim 



log Mg(r) 
logr 



(3) 



For well-behaved multifractals, this function de- 
termines the multifractal spectrum through a Leg- 
endre transform (Falconer, 2003): 



/(a) =infg[ga-r(g)]. 



(4) 



where inf^ means the infimum with respect to q. 
If T{q) is differentiable and convex, the infimum 
is given by the vanishing of the derivative with 
respect to q, so 

"(?) =T'{q)- 
Then, after inverting a{q), one obtains 



f{a) ^ q{a) a - T[q{a)], 



(5) 



which is also convex, namely, f"{a) < 0. 

It is useful to define the quantity D(q) = 
T{q)/{q — 1), which is constant for a monofrac- 
tal, and, in general, decreases with q. It turns out 
that the most important values of q are q = 0, 1. 
At g = 0, we have /(ao) = £1(0) ^ -r(0). On 
account that /'(ao) = 9 = 0, the set of singu- 
larities with exponent ao has the largest fractal 
dimension (which often corresponds to the mea- 
sure's support). The value q = 1 gives the en- 
tropy dimension -D(l), for a value of ai such that 
ai — f{ai) — -D(l). The corresponding singu- 
larity set defines the measure's concentrate, that 
is, the set outside of which the measure vanishes. 
Note that /'(ai) = q{ai) = 1 and convexity of 
f{a) imply /(a) < a for any a (with equality at 
ai). 

For a monofractal, it is sufficient to consider 
the mass-radius relation M2{r) and its exponent 
t(2) — a{q) — D{q) for all q. In finite samples 
from fractals, it is usual to consider the discrete 
version of the mass-radius relation, namely, the 
number-radius relation N{r), given by the number 
of particles in a ball of radius r centered on one 
particle and averaged over every particle. 

In general, the mass-radius relation M2{r) has 
exponent t(2) — D{2). According to Eq. (5) and 
the condition f{a) < a. 



f{a2) = 2a2 



r(2) <a2-^ /(aa) < "2 < D{2), 



where a2 = a{q = 2). We can interpret these rela- 
tions as follows. The multifractal analysis through 
moments relies on a basic property: only concen- 
trations with one exponent contribute to one q- 
moment, namely, concentrations with the expo- 
nent a{q). Therefore, to the mass-radius relation 
only contribute concentrations with exponent a2. 
These concentrations are distributed in a fractal of 
dimension /(a2). The exponent a2 is larger than 
/(a2) due to the presence of other singularities 
with different a in its neighbourhood. Moreover, 
since M2(r) is obtained by averaging over points, 
this operation further increases the exponent. 

2.1.1. Coarse analysis and counts-in- cells 

The mathematical definition of q-moments pro- 
vided by the integral in Eq. (2) may not be conve- 
nient for calculations. Therefore, in analogy with 
the box-counting method for fractal sets, there is a 
method for coarse multifractal analysis (Falconer, 
2003): one puts an ^-mesh of cubes in the total 
volume (which is itself a large cube) and considers 
the mass inside each cube. Then one defines 



MS) 



^^1 



(6) 



where the sum is over non-empty cubes. This 
method is related to the coarse-graining method 
used in cosmology for finite samples under the 
name of "counts-in-cells" (Balian & Schaeffer, 
1989- A): assuming that all the particles have the 
same mass, the mass is measured by counting par- 
ticles. In coarse multifractal analysis, the linear 
cube size £ defines a scale that one lets go to zero, 
eventually. In particular, one lets £ ^ to define 
the function T{q) according to Eq. (3) after replac- 
ing r with L From a mathematical standpoint, the 
equivalence of both point-centered results, from 
Eq. (2), and lattice-based results, from Eq. (6), 
is not guaranteed and requires a careful analysis 
(Falconer, 2003; Hartc, 2001). 

Since every set of points in which a takes a 
definite value is a fractal with dimension f{a), in 
coarse multifractal analysis, the number of cells 
with exponent a is given by 



N{a)^rf^°'\ 



(7) 



scription provided by its number-radius relation 

N{a,r) -r^("), 

namely, the average number of cells with expo- 
nent a in a ball of radius r ^ i centered on one 
cell. Note that the opposite sign of the exponent 
in these two numbers is in accord with their be- 
ing related to Mq and M2, respectively, and the 
equality t(2) = -D(O) = —'''(O) for a monofractal. 

From relation (7) , we can obtain the mass func- 
tion at scale £, namely, the number of cells with 
mass between m and m + dm, which is 

N{m)dm i-f^^Ua. 

Given that a ^ logm/log£, 

N{m) - r-^'(")-". (8) 

Any measurement of positions of objects (galax- 
ies, etc.) has to be made with a given preci- 
sion, which actually defines the natural cell size 
for the corresponding observation. Furthermore, 
a single object at the given resolution may ap- 
pear as a compound object at a higher resolu- 
tion. Therefore, coarse multifractal analysis is 
naturally adapted to this idea of change in the 
description of a distribution in terms of single ob- 
jects as the resolution increases. However, let us 
note that other coarse-graining procedures and, in 
general, other methods can be employed for mul- 
tifractal analysis. In particular, let us point out 
that Halsey et al (1986) define multifractal anal- 
ysis by coarse-graining into arbitrary (disjoint) 
pieces Si, ... , Sn, such that their diameters fulfill 
£i < i. Then, the function 



Mg{T,{S,},£) 



N q 

2—1 



according to the box-counting interpretation of the 
dimension f{a). Let us mention the alternate de- 



for small £ (large N) and optimized over {Si], is 
of order unity only for some value of r, for given 
q. In the case of partitioning into an £-mesh of 
cubes, this definition of T{q) coincides with the 
value given by Eq. (6). 

2.2. Linear multifractal spectrum: bifrac- 
tals 

The Legendre spectrum given by Eq. (4) is only 
an upper bound to the real multifractal spectrum 



(Falconer, 2003). Nevertheless, for "reasonable" 
multifractals they coincide. However, there is an 
interesting situation in which there can be a dis- 
crepancy, as has been noted by Balian and Scha- 
effer (1989-B). In general, the Legendre spectrum 
(4) must be convex (from above), but the real 
multifractal spectrum needs not be convex. So 
the Legendre spectrum associated to a non-convex 
multifractal spectrum is only its convex hull. 

For example, let us consider a bifractal, with 
two exponents ai and a2 (Balian and Schacffcr, 
1989-B). We assume that ai < a2. Although 
the multifractal spectrum seems to consist of just 
two values, namely, /(ai) and /(a2), the Legendre 
spectrum is a full segment with ends at (ai, /(ai)) 
and (q;2, 7(0:2)), that is, a linear spectrum. The 
corresponding slope of T{q) has a jump from ai to 
a2 at 

_ 1(^2) - .f{ai) 
a2 — ai 
Now imagine that we superimpose on this bifractal 
another bifractal, either {ai, 03} or {0:3, q;2}j such 
that the exponent a^ fulfills ai < a^ < 012 and its 
corrsponding dimension fulfills 



/(as) < 



f{a2) -/(ai) 



012 — 0.1 



— (a3-ai) + J{ai)\ 



namely, such that the point (03,/ (03)) is be- 
low the segment with ends at [ai, f{ai)) and 
(02,7(02)). This second bifractal will be invisi- 
ble in both the function T{q) (which gives the mo- 
ments) and the Legendre spectrum. 

We can consider a linear multifractal spectrum, 
with /"(o) = 0, as a limit of general convex spec- 
tra, with /"(o) < 0. To be more precise, if we 
are given two exponents ai and 012, and their re- 
spective dimensions /(oi) and /(02), convexity 
implies that /(o), for ai < o < 02, lies above the 
segment with ends at /(oi) and /(02). Assuming 
that r(g) is differentiable, r'(g) = o is an decreas- 
ing function that interpolates from 02 to Oi in 
a range [gi,g2]- As qi -^ q2, T'{q) becomes dis- 
continuous and T{q) has an angle, at some value 
g* e [gi,(72]- Then /(a) collapses to the hnear 
spectrum (from above). In this limit, the multi- 
fractal becomes dominated by the two fractals at 
the ends of the segment, with dimensions /(oi) 
and /(02). In other words, this is a bifractal limit 
of the multifractal spectrum. 



We can distinguish two cases, according to 
whether g* is larger or smaller than one. In 
the former case, the measure's concentrate fulfills 
t'(1) = a2 > oi, so 02 = /(02); in the latter case, 
the measure's concentrate fulfills t'(1) = oi, so 
oi = /(oi). In words, we distinguish whether the 
measure concentrates in either the less or the more 
singular fractal set. In the former case, the mea- 
sure's concentrate and support coincide, whereas 
in the latter case, the measure's concentrate cor- 
responds to oi and the measure's support to 02- 
Of course, we can also consider the borderline 
case with q^ = I. In this case, a — f{a) for 
oi < o < 02, so the mesasure concentrates uni- 
formly over the full interval. 

Let us mention three interesting examples of 
linear multifractal spectrum, corresponding to the 
three above-mentioned cases. The first example is 
a random distribution of a finite number of power- 
law mass concentrations with exponent a < 3 
{one concentration suffices). This is a bifractal 
with exponents a and 3, and dimensions /(o) = 
and /(3) = 3, respectively (Gaite, 2005-A). Of 
course, the measure concentrates in the regular 
component. The second example is provided by 
the one-dimensional adhesion model. Reasoning 
about the properties of devil's staircases, She, Au- 
rell & Frisch (1992) and Aurell et al (1997) show 
that the mass distribution produced by that model 
is bifractal: it has singularities of delta-function 
type, that is, with a = f{a) = 0, and there ap- 
pears, in addition, a component with a > I and 
/(o) = 1 (homogeneously distributed, in one di- 
mension) . The extension of these results to higher 
dimension is a moot point (Vergassola et al, 1994). 
The third example corresponds to the borderline 
type bifractal: it is a bifractal formed by the union 
of two separate monofractals (with different di- 
mension). In it, the measure is uniformly concen- 
trated, such that each fractal holds a finite fraction 
of the total mass. 

Let us finally note that any piece of a convex 
multifractal spectrum that is not very curved can 
be approximated by a linear spectrum and, there- 
fore, by the corresponding bifractal limit. This 
property will be useful for the connection with the 
Press-Schechter formalism (Sect. 3.2). 



2.3. Halos as mass concentrations 

Halo models of large-scale structure assume 
that the dark matter is in the form of collapsed 
halos with definite density profile. This profile is 
established by averaging over angular directions, 
hence obtaining a radial profile, whose functional 
form is essentially the same for every halo (but 
with different parameters). Several forms have 
been proposed, the simplest one being a power 
law p{r) oc r~^ , which is also a limit of other more 
complicated profiles. A power law halo profile can 
be related with a power law correlation function 
(Peebles, 1974; McClelland & Silk, 1977; Sheth & 
Jain, 1997; Murante et al, 1997). It is also intrin- 
sically associated with self-gravity: the Newton- 
Poisson equations are scale invariant; so we expect 
that the mass concentrations produced by them 
are also scale invariant. This argument is valid in 
the strongly nonlinear regime, where the dynam- 
ics dominates over the initial conditions. In par- 
ticular, the stable clustering model (Peebles, 1980) 
supports this hypothesis (Sheth & Jain, 1997). 

Therefore, we define halos to have a power-law 
profile, namely, the singular profile p{r) ex r^^ 
with /3 = 3 — a > 0. But it is obvious that we 
cannot assign a definite size to a halo with this 
profile. Indeed, in a multifractal that is scale in- 
variant on ever decreasing scales, as the definition 
requires, it is not possible to assign a size or a 
mass to any structure, because that would break 
scale invariance: any structure must be repeated 
in every size, reaching up to the homogeneity scale, 
which is the only existent scale. In consequence, 
mass concentrations are point-like mass singulari- 
ties and they can only be classified by their singu- 
larity strength, given by the exponent a. However, 
a singular power-law profile cannot be physically 
realized. In fact, what is available is a finite sam- 
ple from a multifractal, which is not scale invariant 
on ever decreasing scales. 

Thus, on small scales, there is a change in the 
nature of the matter distribution, related to the 
breaking of scale invariance, which allows us to 
assign halos definite size and masses. This change 
takes place in the large scale distribution of matter 
and, more clearly, in cosmological A^-body simu- 
lations, as we shall explain. Of course, on scales 
smaller than the scale of this change, the power- 
law halo profile is altered. Its shape may turn 



out to fit a different power law, thus becoming 
a combination of two power laws, like in current 
halo profile models (Cooray & Sheth, 2002). If 
we use coarse multifractal analysis, the breaking 
of scale invariance appears below a given cell size, 
which we choose as the size of coarse-grained ha- 
los. Then, the number of (coarse-grained) halos 
is finite. In addition to the full set of halos, with 
their respective masses, the matter distribution is 
defined by the distribution of their positions. We 
will treat the distributions of halos in greater de- 
tail in Sect. 4. 

Let us note that a uniform coarse-graining 
length is convenient but not mandatory: in multi- 
fractal analysis, one can actually choose pieces of 
very different size, as commented in Sect. 2.1.1. In 
particular, halos can be defined to have different 
size and can be given by a different condition. In 
principle, halos can be defined by a condition on 
the density, as in the usual definition. The possi- 
bility of using various definitions has already been 
noted by Valageas, Lacey & Schacffer (2000). We 
choose a uniform coarse-graining length, following 
Vergassola ct al (1994). It is a convenient choice, 
well adapted to A-body simulations. 

2.4. Voids as mass depletions 

We define a void as a region with very low or 
vanishing density, but we need to make precise 
what this condition means. The multifractal for- 
malism is very useful in this regard. In this formal- 
ism, the singular concentrations to which we asso- 
ciate halos have a power law profile p{r) oc r~^, 
with f] — 3 — a > 0. But, if a > 3, there is no 
singularity. The value a — 3 {f3 — 0) is Si refer- 
ence: a mass concentration with this exponent is 
not really a concentration: its density is regular, 
namely, non-divergent and non-vanishing. Expo- 
nents a > 3 correspond to points of vanishing den- 
sity, which we naturally associate with voids. Fur- 
thermore, the total mass in voids vanishes, since 
the measure's concentrate, with /(a) = a, can 
only be in regions with a < 3 (note that, neces- 
sarily, f{a) < 3). 

However, any neighbourhood of points in voids, 
with a > 3, is not totally empty. It is possible to 
have mass around points with vanishing density. 
There may or may not be totally empty voids, ac- 
cording to whether the measure's support occupies 
the full volume or not: if it does not, there are re- 



gions with non-vanishing volume and no mass in 
them (in mathematical terminology, these regions 
are open sets) . If the dimension of the support of 
the measure is smaller than three, there need to be 
totally empty voids, and the multifractal is called 
lacunar. In particular, there may be a scaling set 
of totally empty voids. This is always the case 
in one dimension if the the measure's support has 
dimension smaller than one, but the geometry of 
voids is more complex in higher dimensions (Gaite, 
2005-B; Gaite, 2006). 

To illustrate the nature of multifractal voids 
in a simple case, let us consider again the one- 
dimensional adhesion model. In this model, the 
mass is concentrated in shocks and their locations 
form a dense set (She, Aurell & Frisch, 1992; Ver- 
gassola et al, 1994). The adjective "dense" is un- 
derstood with its mathematical meaning: a set is 
dense in an interval, say, if in any sub-interval, 
however small, there are points of the set; a typ- 
ical example is the set of rational numbers be- 
tween and 1. In the one-dimensional adhesion 
model, while the density at non-shock points van- 
ishes (they have a > 1), these points are sur- 
rounded by infinitesimally close small shocks that 
make the mass in any neighbourhood of them to be 
non- vanishing. Of course, we could make in that 
distribution a totally empty void by removing the 
mass in a sub-interval. 

With the preceding definition of void, we can 
actually consider a multifractal distribution (in 
three dimensions) as the union of (i) halo clusters, 
that is, clusters of mass concentrations with a < 3, 
(ii) voids, that is, mass depletions with a > 3, and 
(iii) the boundary between both regions, formed 
by points with a — 3. However, we must beware 
once more of relying on our geometrical intuition 
when dealing with multifractals: the set of points 
with a = 3 may not be two-dimensional and so 
it may not constitute a regular boundary. Indeed, 
the dimension of this set is f{a = 3), which can 
well be larger than two. We will see in Sect. 5.3 
what this implies for the shape of voids in an ex- 
ample from a cosmological simulation. 

Let us now consider a coarse-grained multi- 
fractal. The coarse-grained density can be de- 
fined to be regular everywhere (with no singu- 
larities) and it only vanishes at a point if the 
mass vanishes in a neighbourhood of it. There- 
fore, the density only vanishes inside totally empty 



voids. Nevertheless, we can still define voids as 
depleted regions comprising coarse-grained expo- 
nents ai — log mi/ log ^ > 3 (a^ -^ oo inside a 
totally empty void). Note that the boundary of 
voids is then a coarse-grained version of the above 
mentioned fractal surface, so it loses small pieces 
as some singularities get smeared. 

Finally, we notice that the presence of voids in 
the galaxy distribution was a principal motivation 
for the bifractal model of cosmic structure. Let 
us quote Balian & Schaeffer (1989-B): "The ocur- 
rence of a bifractal behaviour is due to the coex- 
istence at the same scales of clusters and voids" . 
According to our discussion of the linear multi- 
fractal spectra in Sect. 2.2, we can deduce that a 
bifractal in which the voids define one of the two 
fractal dimensions must be such that the measure 
concentrates in the more singular set and, in ad- 
dition, such that the more regular set has Q!2 > 3 
and is associated with the voids. An interesting 
one-dimensional example is the adhesion model, 
according to Sect. 2.2. However, "coexistence at 
the same scales of clusters and voids" takes place 
in most multifractals, not only in bifractals. 

3. The density probability distribution 
function 

In coarse multifractal analysis, moments are de- 
fined by Eq. (6). There is a related way to define 
statistical moments, namely, by using the density 
probability distribution function, and this defini- 
tion has been appHed to multifractals (Borgani, 
1993; Valageas, 1999). If P(p) is the probability 
of having density p, then, its moments are 
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Since in multifractals the density does not exist at 
singular points, p in this definition must be also 
understood as a coarse-grained density. Besides, 
this definition implies that /io = 1, assuming P{p) 
normalized. In contrast, Mq 7^ 1, in general, and 
it is related with the fractal dimension of the mea- 
sure's support. In coarse multifractal analysis, the 
precise relationship between both definitions of q- 
moments is 

/; (f^ - ln'i\ - ^jPi - /-Sq ^gffl ^ fTiq}-T{0)-3q 



(where the sums are over non-empty cells). For 
example, fJ,i{£) is the average density (within 
the measure's support) and scales with exponent 
— r(0) — 3 = -D(O) — 3 < 0. If the measure's sup- 
port has dimension 3, then the average density 
is £- independent. In this case, it is natural to 
normalize the density field by dividing it by the 
average density, so that ^i = 1 (Valageas, 1999). 
If the measure's support is fractal, then the aver- 
age density diverges as ^ ^ 0. 

Regarding higher moments /ig(^), q > 1, we 
have 

^gffl ^f{q~l)lD{q)-D{0)] 

{pY 

Given that D{q) is a non-increasing function, this 
quotient is finite, in the limit £ ^ 0, only for a 
uniform multifractal, namely, a regular distribu- 
tion or a monofractal (we say that a monofractal 
is uniform because the exponent a is constant). 
For a generic multifractal, that quotient diverges. 

One can derive other relations between mo- 
ments. For example. 
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^(q-l)[D{q)-D{2)] 



(10) 



which diverges for q > 2, except in uniform mul- 
tifractals, namely, monofractals (including regular 
distributions, with a = 3). Thus, this latter re- 
lation is useful in the context of hierarchical mod- 
els: a scale-independent quotient Mq/M2~ im- 
plies D{q) — D{2). In other words, "hierarchical 
correlation functions imply monofractality of the 
distribution" (Borgani, 1993). An analogous rela- 
tion holds for Hg{£) a D{0)^3. 

The probability distribution function (PDF) of 
the coarse-grained density, Pi{p), is directly re- 
lated to the mass function and, therefore, to the 
coarse multifractal spectrum (Sect. 2.1.1): the 
number of cells with density p is the box-counting 
approximation to the fractal dimension of the 
set of mass concentrations with exponent a = 
log m/ log £ = log p/ log i* + 3. The general state- 
ment of this relation is best made in the context of 
probability theory, namely, of the theory of large 
deviations. The relevant concepts in this theory 
appear naturally when we consider the general- 
ization of the central limit theorem to nonlinear 
distributions and, hence, we regard the lognormal 
distribution as a general nonlinear PDF. 



3.1. Connection with the lognormal model 

As a simple extension of the Gaussian density 
field in the linear regime. Coles & Jones (1991) 
proposed that a lognormal field model would be 
suitable for the nonlinear regime, especially, in the 
weakly nonlinear regime. We shall see that a log- 
normal density PDF is natural in connection with 
multifractals. 

From a purely statistical point of view, the 
Gaussian distribution is associated with the cen- 
tral limit theorem and one expects that it appear 
in many situations where the sum of a large num- 
ber of independent random variables is involved. 
A nonlinear analogue of this theorem is provided 
by the multiplication of independent random vari- 
ables, because multiplication of these variables is 
equivalent to addition of their logarithms. Hence, 
the lognormal distribution arises as the limit dis- 
tribution in multiplicative processes. This formu- 
lation of the lognormal distribution corresponds 
to the large- deviation formulation of multifractals 
(Frisch, 1995; Sornette, 2000; Harte, 2001), which 
is particularly adapted to random cascade pro- 
cesses but has a broader scope. We proceed to 
introducing a few relevant notions, leaving aside 
innecessary details. 

We may recall that the central limit theorem 
considers the limit distribution of the sum of n in- 
dependent random variables with finite variances 
(and zero mean). It states that this limit distribu- 
tion is Gaussian, with variance given by the sum 
of the variances. In the particular case of n in- 
dependent identically distributed variables, their 
sum divided by ^/n has a Gaussian distribution 
(with the common variance). However, if we want 
to take into account unlikely events, that is de- 
viations of magnitude larger than ^/n, we must 
consider higher moments than the variance. In 
particular, if we call X the sum of the n variables, 
as n — > oo. 
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In this formula, x is an "intensive" variable, that 
is, finite in the limit n ^ oo, and S{x) is an en- 
tropy function called the Cramer function. This 
function can be obtained from the full (common) 
distribution function of the variables (Frisch, 1995; 
Sornette, 2000). Its second order expansion at its 
maximum (assumed unique) gives the Gaussian 



approximation that constitutes the central limit 
theorem. 

In multiplicative processes, such as random cas- 
cades or fragmentation processes, we also have 
that a good approximation is given by the two 
basic moments (mean and variance) of the loga- 
rithm of the relevant random variable. Therefore, 
the lognormal distribution appears as the basic 
approximation to the distribution of mass concen- 
trations in a multifractal. This has already been 
noted by Jones, Coles & Martinez (1992). The 
large-deviation formalism becomes necessary if we 
need to consider mass concentrations that appear 
rarely, because of their small or large exponent a 
with respect to ao, the value that maximizes /(a) 
(assumed unique). In fragmentation cascades, the 
number n in the hierarchy gives the scale, for ex- 
ample, £n = 2~" for a binary division. The sum 
variable X represents minus the logarithm of the 
mass, such that the "intensive" variable is now the 
exponent a = logTO/log£„. Therefore, 

P[-log2TO-na] -e"^("\ 

and we are to identify S{a) with the multifractal 
spectrum ,f(a); namely, S(a) — f(a) In 2. Thus 
the large-deviation condition is equivalent to fo- 
cusing on mass concentrations with exponents a 
such that their supporting regions have small di- 
mension /(a) compared with its largest value. In 
the "small-deviation" regime, we can use the sec- 
ond order expansion (parabolic approximation) of 
/(a) at ttoj the value that maximizes it (assumed 
unique). This value corresponds to the measure's 
support. The "small-deviation" regime is indeed 
the weakly nonlinear regime, and it gives a lognor- 
mal PDF. 

Note that Jones, Coles & Martinez (1992) chose 
to take a parabolic approximation to /(a) that is 
not its second order expansion but instead makes 
it tangent to the diagonal, like the full f{a). This 
choice does not seem related to the theory of large 
deviations. However, a large-deviation expansion 
around the measure's concentrate (the point of 
tangency to the diagonal) is possible as well. Let 
us consider the mass moment Mi and its associ- 
ated mass PDF 



m P[— log2 m ^ na\ ^ me' 



'^Sia) n[f(a)-a] In 2 



In the limit n ^ oo, this density vanishes unless 
f{a) = a, that is, unless in the measure's con- 



centrate. In other words, there is no mass outside 
the concentrate of the measure (mass) , according 
with its name. Nevertheless, we can consider again 
a "small-deviation" regime given by the second or- 
der expansion around the measure's concentrate, 
such that it gives a lognormal mass PDF. Note 
that this lognormal mass PDF is generally differ- 
ent from the previous lognormal PDF. 

We must mention that the lognormal approx- 
imation to multifractals assumes that the mul- 
tifractal spectrum f{a) is well behaved, namely, 
that it has second derivative. Of course, this con- 
dition of differentiability may not hold in some 
cases. In particular, a bifractal spectrum does not 
have well-defined second derivative at its maxi- 
mum. 

It is well known that the lognormal distribution 
is not determined by its moments /i„, n — 1,2, . . ., 
that is, there are other distributions with the same 
set of moments. In physical terms, the reason for 
this non-uniqueness can be ascribed to the fact 
that high-order moments are dominated by high- 
density regions, so that the full set of moments 
contains insufficient information about very low- 
density regions or voids (Coles & Jones, 1991). 
This was the motivation for considering the addi- 
tional information on voids, namely, the void prob- 
ability function, as a necessary ingredient. In fact, 
the combination of separate information on clus- 
ters and voids, respectively, and the assumption of 
scale invariance led Balian & Schaeffer (1989-B) to 
propose a bifractal model of large scale structure. 
From the preceding discussion about the lognor- 
mal model and its relation with multifractals, we 
can deduce that Coles & Jones's (1991) and Balian 
& Schaeffer's (1989) arguments are related and 
that we do not need to restrict ourselves to bifrac- 
tals. Indeed, we have already pointed out that the 
singular nature of multifractals requires us to use 
the full set of moments Mq{r), for — oo < g < oo. 
In particular, negative values q < Q are typically 
associated with voids. 

We have explained the connection between mul- 
tifractal and the lognormal model of Coles & Jones 
(1991) in regard to the density PDF. Further con- 
nection in regard to n-point correlation functions, 
if possible, seems to require additional assump- 
tions. We note, in particular, that there is no scale 
invariance in the lognormal model as a field the- 
ory, in spite of the arguments presented by Coles & 
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Jones (1991) regarding a possible connection with 
multifractality. 

3.2. Connection with the Press- Schechter 
mass function 

The usual mass function in halo models is the 
one given by the Press-Schechter spherical col- 
lapse formalism (Mo & White, 1996; Sheth & Jain, 
1997). This formalism relies on an initial Gaussian 
PDF. Regarding the relation that we have estab- 
lished in multifractals between the mass function 
of halos and the coarse density PDF (for a < 3), 
we shall explore if there is a connection with the 
Press-Schechter mass function. Let us briefly re- 
call the relevant elements of this formalism. 

Let us consider the coarse density and its PDF. 
For early times or large values of the coarse- 
graining length £, we can consider the density con- 
stant (/? = 1, if the total volume and mass are 
normalized) except for small fluctuations given by 
a Gaussian PDF. This Gaussian PDF has a scale- 
dependent variance cr^(£) <C 1. To obtain the mass 
function of collapsed objects, one needs a criterion 
to decide which regions of size £ will collapse. The 
simplest criterion is that every region in which the 
density exceeds some threshold pc must collapse. 
Given that the PDF is Gaussian, the cumulative 
distribution function is 



P{P > Pc 



1 



/2na 



exp 



2cr2 



1-erf 



Pc - 1 
V2a 



dp 
(11) 



The mass of regions of size £ is m ^ £^. Therefore, 
we can express a'^{£) in Eq. (11) in terms of the 
mass of collapsed objects, provided that this equa- 
tion gives the fraction of regions of mass greater 
than TO that collapse for a given threshold pc- We 
may call this fraction F^ (m) . It is a function of 
the ratio (pc — 1)/ct, so that raising the threshold 
Pc at fixed a is equivalent to diminishing a{£) at 
fixed Pc and making density peaks smaller. a{£) 
must be a decreasing function of £ and, therefore, 
of m. Thus, collapse of large masses is suppressed. 
The mass function is 



N{m) 



dF>{m) 1 
dm m 
da pc — 1 



1 1 



/27r TO, dm cr^ 



■ exp 



1 dF^jm) da 
n da dm 

jpc-ir 

2a2 



where we have introduced the factor m~^ ~ £~^ 
which represents the total number of collapsed ob- 
jects. We assume an initial power-law spectrum 
of Gaussian fluctuations, namely, a'^{£) ex £^^~^ 
(n > —3 is the spectral index in Fourier space). 
Then, the mass function is also a power law, 
but with an exponential decay for large masses; 
namely. 



AT/ N I ^ 
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n/6-3/2 



exp 



n/3-l-l' 
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(13) 
where to,* stands for the large-mass cutoff. 

We must now consider carefully the limita- 
tions of the Press-Schechter formalism. For this, 
me must distinguish the one-dimensional case 
from the others, including the physical three- 
dimensional case. The essential assumptions are 
the validity of an initial Gaussian PDF and of the 
spherical collapse model. A Gaussian PDF is valid 
if £ is large, namely, a'^{£) ^ 1. The spherical 
collapse model amounts to assuming that every 
overdensity above the threshold gives rise to a 
symmetric collapsed object. This is a reasonable 
assumption in one dimension, but it is question- 
able in higher dimensions. In fact, only a small 
percentage of the initial overdensities undergo 
nearly spherical collapse, and they are precisely 
the largest overdensities (Audit, Teyssier & Alimi, 
1997; Lee & Shandarin, 1998). This is particularly 
clear in the Zeldovich approximation. Thus, the 
Press-Schechter formalism can only apply, in three 
dimensions, to the formation of the strongest con- 
centrations (as already noted by Vergassola et al 
(1994) in their comparison of the Press-Schechter 
formalism with the adhesion model). The Press- 
Schechter formalism has been extended to ellip- 
soidal collapse, with the consequent improvement 
(Sheth & Tormen, 1999). However, this exten- 
sion still excludes irregular collapses, which are 
common for weak concentrations. 

In conclusion, the Press-Schechter formalism 
produces a power-law mass function with a large- 
mass cutoff, but it can only be applied to massive 
halos. Considering this formalism as an approx- 
imation to the full nonlinear halo mass function 
based on the linear theory, we would like to ex- 
plore its relationship with a lognormal PDF, which 
is valid in the weakly nonlinear regime. Then the 
question is if relatively strong mass concentrations 
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can in fact have a power-law mass function (ne- 
glecting the exponential cutoff) that is consistent 
with an overall lognormal-like distribution. 

To answer this question, let us mention firstly 
that it is known that the lognormal distribution 
can be mistaken for a power law over a relatively 
large interval (Sornette, 2000, page 80).^ As a 
parabolic approximation to a multifractal spec- 
trum, the lognormal distribution derives from the 
second order expansion around its maximum, rep- 
resenting the measure's support. But the mea- 
sure's concentrate is best represented by the sec- 
ond order (parabolic) approximation around the 
point f{a) — a, as explained in Sect. 3.1. This 
second parabola can actually be very close to the 
first one. In any event, an even simpler approx- 
imation is the first order (linear) approximation 
around the point f{a) = a. 

In fact, a poor sampling of a multifractal will 
mainly reflect the measure's concentrate. Let 
tti = q;(1); given that ai = /(ai) and /'(ai) = 1, 
the linear approximation around the measure's 
concentrate has unit slope, such that it becomes 
an even bifractal (not loaded toward either the 
more or the less singular end, as explained in Sect. 
2.2). In contrast, note that the power-law part of 
the Press-Schechter mass function, Eq. (13), pro- 
duces a divergence in its integral over m in the 
limit TO ^- oo. Of course, this divergence is killed 
by the exponential factor. Nevertheless, it indi- 
cates that the mass concentrates on massive ha- 
los. In this regard, the Press-Schechter mass func- 
tion is analogous to the mass function of the one- 
dimensional adhesion model (according to our dis- 
cussion in Sect. 2.2). This is not surprising: the 
one-dimensional adhesion model is the simplest 
mathematical description of the notion of collapse 
of overdensities, realizing a collapse into delta- type 
singularities (Vergassola et al, 1994). 

4. The fractal distribution of halos 

The realization of scale invariance in a multi- 
fractal is characterized by the multifractal spec- 



^This fact has also been noted in the astrophysical htcra- 
ture, in a small-scale process: star formation from cloud 
fragmentation. The theory of Zinnecker (1984) leads to a 
lognormal mass function, as mentioned by Coles & Jones 
(1991). Zinnecker (1984) compares the lognormal mass 
function with Salpeter's power law. 



trum. In Sect. 2.1.1, we have shown its relation to 
the mass function, which is related in turn to the 
PDF Pe,{p). This PDF can be represented by the 
statistical moments ^q{t} (Sect. 3). However, in 
cosmology, it is normal to measure correlations of 
a given class of objects rather than moments of the 
total dark matter distribution. More precisely, one 
measures positions of luminous objects, namely, 
galaxies, and computes their correlation functions 
to establish the statistical properties of the distri- 
bution, which may include scale invariance, etc. 
One can associate galaxies with dark matter halos 
and, therefore, deduce statistical properties of the 
distribution of these halos; and viceversa. This is 
why the distribution of halos has special interest. 

If the full dark matter distribution is scale in- 
variant, then the correlation functions of the dis- 
tribution of halos must be power laws. More pre- 
cisely, in a multifractal distribution, every popu- 
lation formed by similar halos is a monofractal, 
although different populations have different di- 
mensions (Gaite, 2005- A). We can describe this 
difference between populations as a kind of bias, 
albeit of non-linear type. 

Actually, the very definition of multifractal im- 
plies its interpretation as a set of fractal distri- 
butions of mass concentrations (and depletions), 
classified according to the different values of the 
exponent a: a set of concentrations with expo- 
nent a forms a monofractal of dimension f{a). 
However, for a given exponent a*, the correspond- 
ing set of concentrations generally has dimen- 
sion /(a*) < a*, whereas a monofractal fulfills 
the equality a* — f{a*). This happens because 
in the neighbourhood of the concentrations with 
a = /(a*) are present other concentrations with 
a 7^ /(a*), the mass of which raises the singularity 
exponent from /(a*) to a*. 

In the coarse-grained formulation with fixed 
halo size £h, we have a finite number of individ- 
ual halos which we can classify by their masses. 
Then, each class defines a population and we can 
compute its number-radius relation, N{a, r), a ^ 
logm/logih (Sect. 2.1.1). The fractal dimension 
obtained from each number-radius relation is the 
only quantity necessary to determine the corre- 
sponding monofractal scaling. Thus, the full set of 
dimensions obtained in this way constitutes an ap- 
proximation to the multifractal spectrum (in the 
range a < 3). 
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We have defined halos by relying on coarse niul- 
tifractal analysis with a given cell size. In con- 
trast, the measure N{a,r) is point-centered, like 
the mass-radius relation, given by Eq. (2) with 
q = 2. However, this combination of lattice and 
point-ccntcrcd methods is consistent. Neverthe- 
less, one can also use the coarse-grained probabil- 
ity Pi{p), providing information on the number of 
boxes of size i with mass ?n ~ £" and, hence, on 
the multifractal spectrum (as explained in Sect. 
3.1). Using several values of £, we would repro- 
duce the histogram method of calculating the mul- 
tifractal spectrum (Falconer, 2003), which indeed 
consists of making histograms of the ai for several 
i. This method is, in general, less suitable than 
the method of moments (which uses the Mq{i) of 
Eq. (6)). In fact, ee have found most convenient 
the procedure that consists of (i) the definition of 
halo populations at given ih, and (ii) the calcu- 
lation of their respective number-radius relations 
N{a, r). The explanation lies on the effect on cor- 
reslations of a limited scaling range, as we show 
next. 

4.1. The scaling of correlation moments 

In general, the mass distribution can be de- 
scribed by either the coarse-grained probability 
Pi{p) or the moments Mq{r) (or Pq{t)). In partic- 
ular, its multifractal properties can be described 
by either /(a) or T{q), which are related by the 
Legendre transform (Sect. 2.1) and provide iden- 
tical information, in principle. The formulation 
of a multifractal as a fractal distribution of halos 
emphasizes the scale invariance of each halo pop- 
ulation; in particular, the number-radius relation 
N{a,r) expresses the scaling properties of each 
halo population as a monofractal, such that the 
full set of exponents of all populations constitutes 
an approximation to the multifractal spectrum (in 
the range a < 3). Therefore, the full set of expo- 
nents also provides an approximation to T(g), in 
the corresponding range of q. 

The full mass-radius relation M2(r) is directly 
related with the reduced two point-correlation 
^(r), such that scale invariance implies that they 
both are power laws. Unfortunately, the analysis 
of the two-point correlation function in cosmolog- 
ical iV-body simulations by the Virgo Consortium 
has not shown any convincing scaling (Jenkins et 
al, 1998). We do not consider this negative result 



as ruling out scale invariance. In fact, we can un- 
derstand it as an effect of a limited scaling range in 
a distribution that is actually multifractal, rather 
than monofractal. Let us explain how that effect 
arises. 

The equivalence between f{a) and T{q) relies 
on the one-to-one dependence a{q). However, this 
one-to-one dependence only holds in a mathemat- 
ical multifractal (in the limit of ever decreasing 
scales). In the coarse-grained formulation, we 
can easily deduce that the average over popula- 
tions implied by the definition of Mq{r), Eq. (2), 
will spoil the scaling of each particular population 
(with a given exponent a). In other words, in the 
coarse-grained formulation, every halo population 
contributes to some extent to a g-momcnt (in par- 
ticular, to M2(r)), rather than only the popula- 
tion with a{q). For this reason, to realize scale 
invariance, it is preferable to analyse the fractal 
distributions of halos independently, rather than 
the g-momcnts of the total distribution. 

4.2. The transition to homogeneity 

So far, we have used the definitions of scaling 
exponents, for example, Eq. (1), as impliying naive 
scaling relations, such as to ^ r". Of course, this 
relation holds in a limited range of scales (namely, 
of logr). In addition to the restriction of scaling 
on small scales, already studied, there has to be 
a transition to homogeneity on large scales. Over 
the scale of this transition, say tq, the relation 
m ^ r°' undergoes a crossover to ?Ti ^ r^ (ho- 
mogeneity). Therefore, we can define tq by the 
equation 



A-Kp 3 



ro 



r <^ro . 



(14) 



where p is the large-scale mean density. Analo- 
gously, the halo number-radius relation N(a,r) ^ 



rfi^) 



(Sect. 2.1.1) can also be written as 



N{a,r) 



47rn(a) 



^'o 



/(") 



, r<ro, (15) 



where n{a) is the large-scale mean density of halos 
of exponent a. 

The procedure to determine rg from N{a,r) 
consists of an analysis of the crossover to homo- 
geneity in its log-log graph. It must be done 
for several values of a, to check that the results 
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coincide. Actually, the extreme cases are most 
useful; that is to say, we may take, on the one 
hand, (a) the smallest value of a (the largest 
halo mass) such that there is a sufficient number 
of halos to compute N{a,r), and, on the other 
hand, (b) the largest value of a (the smallest halo 
mass) such that its fractal dimension is sufficiently 
smaller than three as to let us clearly perceive the 
crossover to homogeneity. In Sect. 5.2, we shall see 
how to carry out this procedure in an example. 

5. Halos and voids in cosmological simula- 
tions 

In A'^-body cosmological simulations, the nature 
of the distribution changes on small scales because 
there is (i) a minimal mass, corresponding to one 
particle, and (ii) a softening of the gravitational 
force. The minimal mass defines the discretization 
scale, namely, the linear size of the volume per par- 
ticle. The larger of the discretization scale and the 
softening scale, usually, the discretization scale, is 
the natural (intrinsic) coarse-graining scale of the 
distribution. 

Since all the particles have the same mass, the 
mass is measured by counting particles. We have 
already defined coarse multifractal analysis and 
mentioned its relation with the method of "counts- 
in-cells" (Bahan & Schaeffer, 1989-A), whereby 
one puts a ^-mesh of cubes in the total volume 
(which is itself a large cube) and measures the 
mass inside each cube (cell) by counting particles. 
In this method, initially and as long as the evolu- 
tion is linear, cells of the size of the discretization 
scale contain one particle per cell. Therefore, ha- 
los only arise in the nonlinear stage, as they con- 
centrate particles from other regions that become 
voids. 

We define mass concentrations (or depletions) 
at scale £ and measure their exponents by the 
equation ai — logrm/ log£, where £ is the cell size. 
Standard coarse multifractal analysis proceeds to 
either (i) make histograms of the ai for several i 
(the histogram method) or (ii) calculate the mo- 
ment sums Mqi£) in Eq. (6) (Falconer, 2003). Ei- 
ther way leads to an estimation of the multifrac- 
tal spectrum. However, here we shall combine 
the "counts-in-cells" method with point-centered 
methods, according to Sect. 4. 

We have applied our numerical methods to sev- 



eral A^-body cosmological simulations. We present 
the results of the analysis of the redshift z — 
positions of the ACDM GIF2 simulation (by the 
Virgo Consortium), with 400^ particles in a vol- 
ume of (110 h~^ Mpc)^. The force-softening 
length is e — 7 h^^ kpc. This simulation is de- 
scribed by Gao et al (2004), who also use it for 
halo analyses, though they employ the usual FoF 
definition. Naturally, the discretization length is 
400~^ (in box- length units) and is considerably 
larger than the softening length. We determine 
from it the appropriate halo coarse-graining scale: 
since we operate with powers of two, we take as 
halo size £h = 256"-'^. Actually, our results are not 
very sensitive to the precise value of £h- 

5.1. Counts in cells and halos 

In the method of "counts-in-cells" , the qualita- 
tive effect of discreteness dependes on the cell size, 
as we study next. If this size is too small, namely, 
much smaller than the force-softening length, most 
cells are empty an the non-empty ones have only 
one particle, so this is the regime totally domi- 
nated by the discretization and the gravitational 
softening. As the size grows, the ratio of non- 
empty cells grows as well, and some of them begin 
to have more than one particle. In fact, as the cell 
size grows, soon a few cells become massive (while 
the great majority are still empty), in such way 
that a pattern appears: the number of cells with 
given mass N{m) becomes a power law. This is 
shown in Fig. 1.^ This regime is still dominated 
by the discretization. 

As the cell size grows larger, at some point, the 
most numerous cells are no longer the ones with 
TO = 1 (Fig. 1). Further increase of the coarse- 
graining scale eventually leads us to a lognormal- 
like distribution. We note that this last transi- 
tion, namely, from the power-law regime to the 
lognormal-like regime, takes place at a coarse- 
graining length aproximately equal to £h — 256~^ 
(Fig. 1). Thus, this is the scale at which we have 



^Thc power-law dependence of N(m) implies, in particular, 
that the cumulative cell number Ny (m) (the number of 
cells with mass larger than m) also follows a power law. 
Since 7V>(m) is the rank, we have that the rank-ordering 
of cells follows a power law, that is, it satisfies Zipf's law 
(Zipf, 1949; Sornette, 2000). This property also applies to 
halos, according to our definition. The connection of the 
Zipf law for voids with fractality has been studied by Gaite 
& Manrubia (2002) and Gaite (2005-B, 2006). 
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maximal richness and variety of halos, that is, they 
span the maximal mass-scale range. On the con- 
trary, at this scale, voids are hardly sampled. 

Within the power-law regime of N{m), the 
power law holds very clearly. For example, in 
the top plot in Fig. 1, in the horizontal axis (log 
m), the scaling holds from to 7, that is, from 
TO = 1 to m = e^ ~ 1097, namely, three decades 
(and larger range on the vertical axis). We have 
performed least-squares fits of log N versus log m 
to find the power-law exponents. To minimize 
the standard error of the fits, it is convenient 
to reduce the scaling range. The fit correspond- 
ing to £ = 256~^ in the range m e [10, 100] 
yields —1.996 ± 0.004; the fit corresponding to 
£ ^ 512-1 yigijg -2.022 ± 0.004. We conclude 
that N{'m) ~ m~^. This power law amounts to 
a linear multifractal spectrum, according to Eq. 
(8). Moreover, it is consistent with a — ai, corre- 
sponding to the measure's concentrate (such that 
ai = f{o-i) and /'(ai) = 1). 

The counts-in-cells method of computing N{m) 
for various coarse-graining lengths would allow us 
to estimate the multifractal spectrum (according 
to the histogram method). However, the counts- 
in-cells formalism is especially suited to compute 
the moment sums, which provide a more accurate 
estimation. In fact, it is best to replace in Eq. (6) 
for (7-moments the sum over cells with a sum over 
number of particles: 
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log m 



Fig. 1. — Log- log plots of number of halos N ver- 
sus their mass m (number of particles) at coarse- 
graining scales 512"^ (top), 256-^ (middle), and 
128^1 (bottom). Both the top and middle plots 
are linear, but the latter, corresponding to lower 
resolution, already shows a slight decay in the 
number of light halos. The bottom plot clearly 
shows a maximum for to, = 2. 



Mq{i) = J2 N{m)m'i 



(16) 



This formula is a discrete version of Eq. (9) . In a 
cosmological simulation, the upper limit to to, in 
the sum is, of course, the total number of parti- 
cles iVtot (400^ in the GIF2 simulation). In addi- 
tion, we have to introduce a normalization factor 
to have Mi ~ 1; namely, we have to divide m by 

iVtot. 

We can deduce from Eq. (16) how halos of dif- 
ferent mass contribute to q-moments. Since we 
have found that the halo mass function N{m) ^ 
TO"^, g-moment are given by sums of g— 2-powers. 
In these sums, the contribution of each term in- 
creases or decreases with to, according to whether 
q is larger or smaller than 2. In particular, every 
term of M2{i) contributes equally. This blatantly 
contradicts the assumption that only contribute 
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to M2{£) singularities with a given exponent (q!2) 
and, therefore, a given m. 

The assumption that only one a contributes, 
namely, a{q), is crucial for the scaling of moments, 
as commented in Sect. 4.1. We display in Fig. 2 
the log-log plot of M2(£), for -9 < log2^ < -1, 
to show the absence of a convincing scaling range. 
Indeed, to reach full multifractal scaling, we need 
to be in the lognormal-like regime of the halo mass 
function iV(m). But this regime appears for cell 
sizes larger than £h and up to the transition to ho- 
mogeneity, comprising a scale range that is insuffi- 
cient for the scaling of (/-moments. Therefore, it is 
natural that no scaling range of the two-point cor- 
relation function was found in the analyses of iV- 
body simulations by the Virgo Consortium (Jenk- 
ins et al, 1998). In contrast, scaling is found in 
the distributions of halos, as we show next. 

5.2. Fractal distributions of halos 

The description of multifractal large-scale 
structure in terms of fractal distributions of halos 
has been explained in Sect. 4. In brief, given the 
halo size, halos are classified by their masses (or 
by their exponents, which is equivalent). Then we 
can directly test if each class of halos is a fractal 
set by calculating its two-point correlation func- 
tion, or, rather, its number-radius relation. 

The study of fractal distributions of equal-mass 
halos serves two purposes: first, to provide us with 
a proof of multifractality and, second, to obtain 
the homogeneity scale. According to Sect. 4.2, to 
obtain the homogeneity scale, it is convenient to 
select two halo populations with extreme values 
of the halo mass (but the light halo population 
must not be too homogeneous). Actually, we need 
a range of masses in each population to have a 
sufficient number of halos to compute the number 
radius relation. We choose two ranges, with 750 
to 1000 particles and with 100 to 150 particles, 
defining two populations consisting of 2508 and 
25556 members, respectively (see Fig. 1, middle 
plot). 

Results of the analysis appear in Fig. 3: on the 
top, there is a 1/8-thick slice that shows the aspect 
of both halo populations; below we have log-log 
plots of the number-radius functions for the re- 
spective spatial distributions. They have similar 
scaling ranges, but corresponding to quite differ- 



ent fractal dimensions, namely, 1.1 and 1.9. The 
scale of transition to homogeneity, according to 
Eq. (15), is tq ~ 1/4. Intermediate mass values 
(^ 300) yield intermediate dimensions but similar 
ro- 

5.3. The structure of voids 

In our definition of void, in Sect. 2.4, we out- 
lined how coarse-graining affects the structure of 
voids. As said above, in A^-body simulations, ha- 
los only arise in the nonlinear stage, as they con- 
centrate particles from other regions that become 
voids. Moreover, these voids are totally empty, 
because a cell with less than one particle has to 
be empty. However, an improvement in mass res- 
olution can be achieved with a relatively simple 
numerical technique: re-simulation of selected re- 
gions. Indeed, re-simulation of voids with higher 
resolution unveils new structure, similar to the 
structure with the initial resolution (Gottlober et 
al, 2003). This new structure includes small ha- 
los. In the words of Gottlober et al (2003), "the 
haloes are arranged in a pattern, which looks like 
a miniature Universe" . This type of self-similar 
structure corresponds to a multifractal. 

In our analysis of the ACDM GIF2 simulation, 
we have taken £h = 256, so the smallest halos 
have about 4 particles and cells with fewer parti- 
cles belong to voids. We note that the boundary 
of voids cannot have a good pictorial definition in 
this situation. To appreciate the complicated ge- 
ometry of this boundary, we have plotted in Fig. 4 
a smoothed version of a two-dimensional section 
of it (obtained by means of isodensity contours). 

The voids displayed in Fig. 4 look quite empty. 
However, some particles lie in them. Moreover, 
re-simulation of a region in them with higher res- 
olution (as in Gottlober et al, 2003) would surely 
unveil new structure similar to the structure al- 
ready formed by halo clusters in Fig. 4, but in 
smaller size. That is to say, many small halo clus- 
ters would pop out in the voids. At the same time, 
if the re-simulation were extended to the roundish 
small halo clusters already displayed in Fig. 4, they 
would become less round and more crumpled, like 
the large ones. 
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5.4. Multifractal spectrum 

Finding fractal distributions of halos with defi- 
nite scaling exponents constitutes a proof of mul- 
tifractality. However, we must also apply the stan- 
dard method of multifractal analysis, namely, the 
calculation of g-moments and, hence, the multi- 
fractal spectrum, according to Eqs. (3) and (5). 
Since g- moments have no well defined scaling, T{q) 
is not well defined. Nevertheless, it can be com- 
puted at any given scale, understanding Eq. (3) 
as a an equation at the given scale rather than a 
limit. Then, we can study how the multifractal 
spectrum given by Eqs. (4) and (5) changes across 
the scales. 

Therefore, we have computed the multifractal 
spectrum from g-moments for various £„ = 2~", 
namely, logj £ = —9, . . . , —5, and plot the results 
in Fig. 5. We can see that the multifractal spec- 
trum is stable in the range of a where all the 
curves overlap. For the decreasing part of f{a), 
which corresponds to negative q, we have only two 
curves. Their agreement is sufficient for smaller a 
but they separate for larger a (a > 4). This was 
to be expected, for these values correspond to the 
most depleted voids, which are not sampled. 

The collapse of the f{a) curves along such a 
considerable range of a and for a considerable 
range of scales is surprising indeed (compare with 
the log-log plot of M2{i) in Fig. 2). Furthermore, 
the stability of f{a) contrasts with the change of 
the halo mass function from power law to lognor- 
mal that takes place over the smaller scales in Fig. 
5. In this regard, our conclusion is that the cru- 
cial condition in the calculation of the multifractal 
spectrum by the method of moments has some sort 
of self-consistency; namely, the condition that only 
singularities with exponent a{q) contribute to the 
g-moment (Sect. 4.1), allows one to extrapolate the 
multifractal spectrum to small scales. 

Here is a number of consequences that can 
be directly extracted from the coarse multifractal 
spectrum, such as is displayed in Fig. 5: 

• The entropy dimension is ~ 2.5. This means 
that this multifractal is concentrated in a set 
which is not very far from being homoge- 
neous (dimension three). 

• The largest fractal dimension is actually 
three, so this is the dimension of the mea- 



sure's support. In conjunction with statisti- 
cal homogeneity, this essentially means that 
there are no totally empty voids in the full 
dark matter distribution (see the comments 
in Sect. 5.3). 

• a{q = 0) ~ 3.3 > 3. This value corre- 
sponds again to the multifractal support. 
We remark that it is formed by mass deple- 
tions, that is, such that they belong to voids. 
Therefore, voids (though not totally empty) 
dominate the distribution (see also Fig. 4). 

• The set of regular points with a — 3 (non- 
vanishing density), namely, the boundary of 
voids, turns out to have f{a) ~ 2.9. This 
fractal dimension is very close to three. This 
does not necessarily mean that each piece of 
surface is very crumpled, because the fractal 
dimension is also a measure of fragmentation 
(Mandelbrot, 1977). Therefore, its large di- 
mension surely results from the proliferation 
of small halo clusters that fill the voids. This 
type of geometry can be inferred from Fig. 4. 

6. Conclusions 

We have shown that a statistically self-similar 
multifractal model explains the relevant features 
of the geometry of the large scale structure of mat- 
ter. We have described this multifractal model 
as consisting of singular mass concentrations and 
regular mass depletions (plus a boundary between 
both). Hence, we have proposed to identify mass 
concentrations with halos and to identify mass de- 
pletions as belonging to voids. The identification 
of halos with mass concentrations implies, in the 
scale-invariant limit, that they have an average 
power-law profile, although they can actually be 
very anisotropic. In addition, halos are clustered 
with a power-law form of its two-point correlation 
functions. In particular, every population of mass 
concentrations with equal strength (singularity ex- 
ponent) has a definite fractal dimension. 

We have emphasized that scale invariance, in 
fact, forbids one to distinguish the mass distribu- 
tions of matter inside halos from the distribution 
outside. That is to say, scale invariance actually 
deprives halos of meaning as entities with a given 
size (or mass), so they can only be classified by 
their strength exponent. Therefore, halos can only 
be given a size in connection with the breakdown 
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of scale invariance. In this connection, the number 
of halos becomes finite and their average profile 
changes. In cosmological TV-body simulations, we 
have argued that it is most convenient to use the 
uniform coarse-graining length given by the dis- 
cretization scale (the linear size of the volume per 
particle). 

The uniform size of halos in our multifrac- 
tal model can be interpreted in terms of stan- 
dard coarse multifractal analysis. In TV-body sim- 
ulations, which use equal-mass particles, coarse 
multifractal analysis becomes equivalent to the 
method of counts-in-cells. Each halo is formed by 
a variable number of particles, a number which can 
be small. Each population of halos with (about) 
the same number of particles forms a monofractal 
with a definite dimension. We have analysed, in 
the Virgo A-CDM GIF2 simulation, the mass func- 
tion of halos and their distribution. The distribu- 
tions of the different populations of halos show 
well-defined scaling ranges, corresponding to dif- 
ferent scaling dimensions, and shows a common 
scale of transition to homogeneity. We have se- 
lected two representative populations of heavy or 
light halos, resulting in fractal dimensions D = 1.1 
or D = 1.9, respectively. Note that these two 
values span the usual range of fractal dimensions 
measured in the galaxy distribution. The conclu- 
sion of our analysis of the GIF2 simulation is that 
scale invariance holds in the nonlinear regime re- 
produced in iV-body simulations and that it can 
be measured, contrary to the conclusion of Jenkins 
et al (1998). 

The scale of transition to homogeneity in the 
fractal distributions of halos of the GIF2 simula- 
tion is, in physical units, tq ~ 14ft,~^ Mpc. Al- 
though this value corresponds to a dark-matter 
simulation, it roughly agrees with standard values 
found in the galaxy distribution. 

The fractal distributions of halos with differ- 
ent dimensions are biased with respect to the full 
particle distribution, which has fixed mass-radius 
dimension D{2). Of course, this bias is different 
from the usual linear bias assumed in cosmology, 
which only affects the amplitude of the two-point 
correlation function: it is a bias that affects the ex- 
ponent, that is to say, the fractal dimension. We 
may properly call this kind of bias nonlinear bias 
of scaling type. 

Mass depletions, belonging to voids, are very 



undersampled in A^-body simulations: by defini- 
tion, depleted halo-size cells have no particles; or 
they have a few of them, at the most, if we de- 
fine the halo size to be somewhat larger than the 
discretization scale. So the few particles in voids 
are a gross representation of the small halos that 
one would find there at a higher resolution. In 
this regard, our conclusion is similar to that of 
Kuhlman, Melott & Shandarin (1996). However, 
self-similarity allows us to draw further conclu- 
sions regarding the distribution of matter in voids, 
since the properties of the distribution can be ob- 
tained by coarse multifractal analysis with larger 
cell sizes. 

As to the role of correlation moments, multi- 
fractals are defined by the full set of moments 
Mq{r) for — oo < q < oo (or the PDF moments 
Hq{r) plus Mo(r)), whereas mo no fractals are de- 
fined only by M2{r). In particular, monofrac- 



tals fulfill the hierarchical relation M„ 
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for any q. Thus, the singular and heterogeneous 
nature of multifractals implies that the moments 
with q — 2,3,... are not sufficient. In fact, most 
important are q — 0,1, associated to the the 
measure's support and concentrate, respectively. 
Larger g's are also relevant, in particular, q = 2, 
but the information gathered from large-g mo- 
ments only applies to the strongest concentrations 
(most massive halos). The ncgative-g moments 
are related to the distribution of matter in voids. 
Therefore, they are affected by the undersampling 
of voids. However, if the scaling range is suffi- 
ciently large, these negative-g moments can be cal- 
culated (approximately), as well as the part of the 
multifractal spectrum given by them. 

Regarding (/-moments of the full mass distri- 
bution in A^-body simulations, we have seen that 
they scale poorly. In contrast, the analysis in 
terms of fractal populations of equal- mass halos by 
their number-radius relation (or two-point correla- 
tion function of their positions) provides good scal- 
ing, suitable for calculating the fractal dimension 
of each population. The reason for the poor scal- 
ing of (/-moments is that the one-to-one correspon- 
dence between mass concentrations (or depletions) 
and g- moments given by the function a{q) is only 
approximate when the scaling range is limited. If 
several concentration strengths contribute to the 
same g-moment, their mixing spoils the scaling of 
that (/-moment. We have discussed in Sect. 5.1 
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how this mixing of halo populations takes place. 

In particular, the various halo populations con- 
tribute evenly to M2 and, hence, to the full two- 
point correlation, altering its scaling. Therefore, it 
is natural that Jenkins et al (1998) observed devia- 
tions from scaling in their analysis of A^-body sim- 
ulations. It is also natural that Zehavi et al (2004) 
observe deviations from scaling in their analysis 
of the Sloan Digital Sky Survey: if one associates 
galaxies to halos, their distribution is likely mul- 
tifractal, so the mixing of populations takes place 
as well. In fact, Zehavi et al (2005) and Tikhonov 
(2006) have found that the slope 7 of the log-log 
plot of the two-point correlation decreases with 
luminosity, which agrees with a multifractal dis- 
tribution. 

In the GIF2 simulation, at the scale of the halo 
size f-h, the mass function is well described by a 
power law, so that it corresponds to a linear mul- 
tifractal spectrum and the distribution is actually 
bifractal. To be precise, we find N(ra) ~ m~'^, 
very accurately. Therefore, the corresponding 
bifractal fulfills a = f{a) and is an even bifractal, 
in which the measure is uniformly concentrated 
(not loaded toward either the more or the less sin- 
gular end, as explained in Sect. 2.2). 

Thus, the bifractal mass function furnished by 
equal-size halos must be distinguished from the 
Press-Schechter mass function, which (apart for 
the large-mass cutoff) is equivalent to a linear 
mutifractal spectrum concentrated on the singular 
end, that is to say, on the heavy halos. Note that 
Eq. (13) shows that a mass function N{m) ^ to~^ 
is associated with the spectral index n = —3, 
which is just beyond the allowed range. Of course, 
equal-size and Press-Schechter halo mass functions 
were expected to differ. However, it is remarkable 
that both have a power-law range and that the 
constant exponent that appears with our defini- 
tion is a limit case of the Press-Schechter mass 
function. In any event, we must point out that 
the exponent —2 for equal-size halos is indepen- 
dent of the initial power spectrum, in accord with 
the nonlinear nature of the multifractal. 

The mass function of halos is bifractal (a power 
law) , but the calculation of the multifractal spec- 
trum through moments for cell sizes larger than th 
shows convergence towards a convex, lognormal- 
like spectrum (Sect. 5.4). We conclude that, at the 
halo-size scale iu, the mass discretization does not 



allow for a good sampling of halos and allows for 
almost no sampling of voids. Poorly sampled halos 
mainly reflect the measure's concentrate. There- 
fore, the power law is merely a linear approxima- 
tion to the full mass function near the measure's 
concentrate. In this regard, let us remark that, 
in a multifractal that is scale invariant on ever 
decreasing scales, the total mass concentrates in 
the set such that a — ,f{a). Paradoxically, mass 
concentrations with larger a are too few to hold 
any mass, even all together. On the other hand, 
mass concentrations with smaller a (or mass de- 
pletions) are abundant, but their total mass van- 
ishes nonetheless. 

I am grateful to Liang Gao for kindly supplying 
me with the GIF2 data. My work is supported 
by the "Ramon y Cajal" program and by grant 
BFM2002-01014 of the Ministerio de Educacion y 
Ciencia. 
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Fig. 2. — Log-log plot of Al2{i) (logarithms are to 
base 2 and the total size is normalized to unity). 
Note the reduced, hardly noticeable scaling range 
(the straight line corresponds to homogeneity). 
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Fig. 3. — Selection of two halo populations from 
the GIF2 A^-body simulation: heavy haloes with 
750 to 1000 particles and light haloes with 100 to 
150 particles. (Top) 1/8-thick slice showing heavy 
haloes in red and light haloes in blue. (Down) 
Number-radius relation for each halo population, 
showing fractal dimensions 1.1 and 1.9, respec- 
tively, and a transition to homogeneity in both 
(logarithms are to base 2 and the total size is nor- 
malized to unity). 
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Fig. 4. — A section of the isodensity contour of 
the smoothed density field at the density corre- 
sponding to a = 3. The dark area corresponds 
to mass depletions with a = 3 and the light area 
corresponds to mass concentrations with a = 3, 
namely, to halos. The boundary is a very complex 
curve, with dimension larger than one. 




Fig. 5. — Full multifractal spectrum /(a) for 
log2 i — —9, —8, —7, —6, —5 (yellow, red, green, 
blue, magenta). We have also plotted the diag- 
onal line to show how the multifractal spectrum 
is placed under it and touches it at the measure's 
concentrate dimension a ~ 2.5. 
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